The universe is accelerating. Do we need a new mass scale? 
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We try to address quantitatively the question whether a new mass is needed to fit current su- 
pernovae data. For this purpose, we consider an infra-red modification of gravity that does not 
contain any new mass scale but systematic subleading corrections proportional to the curvature. 
The modifications are of the same type as the one recently derived by enforcing the "Ultra Strong 
Equivalence Principle" (USEP) upon a Friedmann-Lemaitre- Robertson- Walker universe in the pres- 
ence of a scalar field. The distance between two comoving observers is altered by these corrections 
and the observations at high redshift affected at any time during the cosmic evolution. While the 
specific values of the parameters predicted by USEP are ruled out, there are regions of parameter 
space that fit Snla data very well. This allows an interesting possibility to explain the apparent 
cosmic acceleration today without introducing either a dark energy component or a new mass scale. 



I. INTRODUCTION 

During the last decade, several observational probes 
[IHH have confirmed that our universe is undergoing a 
phase of accelerated expansion. Beyond the details of 
specific models, one of the most remarkable aspects of 
such a discovery is the seemingly unavoidable presence 
of a new tiny mass scale in the theory that describes our 
world. 

In the framework of General Relativity (GR), a nega- 
tive pressure component ("dark energy" can account 
for the cosmic acceleration. During the cosmological 
expansion, such a component has to become dominant 
when the average energy density p{t) drops to about its 
present value po (t is the proper time). Thus, dark en- 
ergy Lagrangians typically contain a mass parameter of 

the order of M ~ p^ 4 ~ 10~ 3 eV, that triggers the epoch 
when dark energy starts to dominate 1 . Models of mas- 
sive/modified gravity highlight the problem from a dif- 
ferent perspective. If the graviton is effectively massive, 
the modified dynamics of gravity at large distances can 
provide a mechanism for "self-acceleration" Q and/or of 
filtering for the cosmological constant's zero mode Q. In 
that case the new mass brought into the theory is the 
mass of the graviton 2 , m g , typically of order the Hubble 
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1 The mass parameter can be higher e.g. in quintessence mod- 
els with power-law potentials [jj @] , but at the price of giving 
a milder equation of state which is now severely challenged by 
observations. 

2 Interestingly, and as opposed to, e.g., the mass of scalar 
quintessence fields, such a mass might be protected against — 
and actually made smaller by - radiative corrections [jj. 



constant H , i.e. m g ~ H ~ (po/M^) 1 / 2 ~ lCT 33 eV. 
In f(R) theories of gravity, where the Lagrangian den- 
sity / is a function of the Ricci scalar R, the late-time 
acceleration is realized, again, by introducing a curvature 
scale R c of the order of Hq [Io[ (see also Refs. [HI)- Note 
that, although based on different mechanisms, all mod- 
els follow an analogous pattern 3 : there is a tuned scale 
"hidden" in the theory which becomes effective, "by co- 
incidence" , when the appropriate cosmological quantity 
(the average density p(t) or the Hubble parameter H(t)) 
drops to about its value. 

But is the acceleration of the universe intrinsically im- 
plying the presence of a new mass scale? If we allow the 
possibility of departures from GR at large distances there 
is a logical alternative. High redshift observations have 
the unique property of relating objects (e.g., the observer 
and the supernova) that are placed from each other at a 
relative distance of the order of the average inverse cur- 
vature (roughly, the Hubble length ~ Hq 1 )- Therefore, 
modifying GR in the infrared (IR) at a length scale set by 
the curvature - rather than fixed a priori by a parameter 
- will systematically affect any cosmological observation 
at high redshift, regardless of when such an observation 
takes place and without the need of any external mass 
scale. In other words, we might not need a new mass scale 
because we already have (a dynamical) one, the Hubble 
parameter H (t); the only "coincidence" that we might be 
experiencing is that of observing objects that are placed 
from us as far as the Hubble radius 4 . 



3 The few alternatives to this common pattern include the pro- 
posal that our universe is not homogeneous on large scales [T3 | 
and attempts based on possible non-trivial effects of smaller in- 
homogeneities on the cosmic evolution |13| . 

4 The same circumstance does not apply, for instance, to obser- 
vations within the solar system: typical solar system distances 
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The point of view sketched above is somewhat com- 
pelling, it addresses directly the fine tuning and coinci- 
dence problems, but seems to require a serious revision of 
the current low energy framework for gravity. Any grav- 
itational operator that becomes effective in the infrared, 
on dimensional grounds, has to bring in the Lagrangian 
a mass parameter. Moreover, GR itself is already a ge- 
ometrical deformation of flat space at distances of the 
order of the curvature. What seems to be required is a 
further curvature-dependent subleading effect that sys- 
tematically modifies the geometrical description of GR 
at large distances. 

Recently, a proposal along those lines has been made 
by one of the present authors. The modification upon 
the standard framework is forced by imposing an "Ul- 
tra Stron g" v ersion of the Equivalence Principle (USEP, 
see Refs. |l4l. Il5| for more details). USEP suggests that 
the usual geometric description of spacetime as a metric 
Riemannian manifold might hold only approximately, at 
small distances. Such a conjectured "IR-completion" of 
gravity, in its full generality, represents a major theoret- 
ical challenge. However, it can be tentatively explored 
with a Taylor expansion around GR, by applying USEP 
to a specific GR solution 5 (see Appendix [Af . For a spec- 
tator scalar field in a spatially flat Friedmann-Lcmaitre- 
Robertson- Walker (FLRW) universe, the first-order cor- 
rection to GR is calculable, and few cosmological conse- 
quences are derivable [l5j . 

Consider, as the zeroth-order (GR-) approximation, a 
homogeneous, spatially flat FLRW universe with scale 
factor a(t). It is known that in such a solution the phys- 
ical distance d(t) between a pair of comoving observers 
grows proportionally to a(t); otherwise stated, the co- 
moving distance A = d(t)/a(t) is a constant. The first- 
order correction found in Ref. (l5j modifies such an ex- 
pansion law by a subleading distance-dependent contri- 
bution. As a result, the distance d(t) between two comov- 
ing observers grows as a(t) only when small compared to 
the Hubble length H^ 1 but gets relevant corrections oth- 
erwise. Thus, the scale factor a(t) defines the expansion 
everywhere but only in the local limit, and effectively 
detaches from the expansion on the largest scales. 

The corrected expansion is most easily seen in terms 
of the above defined comoving distance A, which is con- 
stant only in the small distance limit. Its derivative with 
respect to observers' proper time t reads (i~5j 



A = A 



3 (H 2 a 2 y 



higher orders, 



(1) 



are always extremely small with respect to, say, the average 
Weyl curvature. An order of magnitude estimate indeed gives 
(Weyl curvature)" 1 / 2 Rj r (r/3 km) 1 / 2 , where r is the distance 
from the Sun. 

5 In a similar fashion, someone who does not know GR can try to 
expand around some point in Riemann normal coordinates and 
find, in some specific cases, the first GR corrections to flat space. 



and is clearly negligible on sub-Hubble scales. For com- 
pleteness, a basic derivation of Eq. (JlJ is sketched in Ap- 
pendix |XJ The comoving trajectory r(t) of a light ray 
also receives corrections because the modified global ex- 
pansion {1} has to be considered on top of the usual con- 
tribution dr = At /a. In a matter-dominated universe, 
where the Hubble parameter at the redshift z = 1/a — 1 
is given by H(z) = Hq(1 + z) 3 / 2 , we have 



d(H Q r) 



1 



(H rf 



(1 + z) 3 / 2 



(2) 



which has to be solved with initial conditions ?*(0) = 0. 

The above modification, including the factor of 1/4, is 
forced by requiring that USEP applies for a scalar field 
in a spatially fiat FLRW universe [Hj]. The correction 
increases the luminosity distance 



d L {z) = (l + z)r(z), 



(3) 



and therefore it effectively goes in the direction of a uni- 
verse with positive acceleration. However, as the present 
analysis shows as a by-product, the correction given in 
the last term of Eq. ([2]) is too small (of too high order in 
the redshift) to explain Snla data. 

Equation @ is an example of a IR-gcomctrical defor- 
mation that does not contain any mass scale but only 
-ff-dependent subleading terms. In this paper, by study- 
ing a generalized version of @, we attempt to address 
quantitatively the question whether Snla data can be ex- 
plained without introducing any new mass scale. We in- 
clude terms of lower order in the redshift that are needed 
to efficiently reproduce Snla observations. Equation ([T]) 
is generalized as follows: 



A = A 1 XH + A 2 X 2 H 2 a + . . . 

+B 1 \ 2 {Hay + B 2 X 3 (H 2 a 2 y 



(4) 



Note that the above structure of corrections includes ([T]) 
as a special case. In practice, the terms in the above 
expansion rearrange when we calculate the luminosity 
distance. Therefore, for a matter-dominated universe, a 
quite general structure of subleading terms is given by 



d(H r) 



1 



F [H r(l + z) 



1/2 



(5) 



dz (I + zfl 2 

where F(x) is a generic function with F(0) = I: 

F(x) = 1 + ax + (ix 2 +'fx 3 + ... . (6) 

By comparison with (0| we have a = —Ai,/3 = Bi/2 — 
A<i , 7 = B-2 ■ Note that Eq. ([2]) corresponds to a = (3 = 
and 7 — 1/4, while in GR all coefficients are set to zero. 

II. EFFECTIVE DESCRIPTION 

It is possible to obtain analytic solutions of §5§ in some 
restricted cases (e.g., 7 = 0, sec Appendix [Bjl . However, 



3 



it is perhaps more useful to study the effective behavior of 
([5]) at low z. In order to make an easy comparison with 
known parametcrizations of dark energy, we can easily 
express our first two parameters, a and /3, in terms of an 
effective density parameter and a constant effective 
equation of state w s of dark energy, in the presence of 
non-relativistic matter 6 . Such effective parameters fl(| 
are thus defined by 



r(z) 



o H cff (x] 



dx . 



(7) 



where 



h cS (z) = H ^J (i - ngf.)(i + zf + n^.(i + z) 3 (i+ w °«). 

(8) 

By expanding Eq. ([7} at small redshift we find 

H r = z - ^- (1 + OrjVeff) 
z 3 

+ — [5 + 8w eff f2ff E + 3w 2 cS n BE (3n UE - 2)] 

+ .... (9) 

On the other hand, the solution of ([5]) can be expanded 
as 

3z 2 / 2 
H a r = z - — I 1 - -a 

z 3 

+— (15 - 14a + 4a 2 + 8/3) + ... . (10) 

By comparing © and (fTU|) we can relate the two sets 
of parameters: 



B = -Weff [1 + 6w efI (fif& - 1)1 nf E . (11) 

o 

For the ACDM model (w cS = -1) with fij^ = 0.7 it 
follows that a = 1.05, B = —0.74. Of course the above 
expansions are valid only for z <C 1, so it is expected 
that the likelihood analysis including high-rcdshift data 
can give different constraints on the model parameters 
(as we will see later). 

The most important correction that leads to a larger 
comoving distance relative to the Einstein dc Sitter uni- 
verse originates from the a term in Eqs. (JSJ)- (JSJ) , i-e., the 
term A\\H in Eq. Note that, for z -C 1, the correc- 
tion 7 becomes important only for the terms higher than 
order z 3 in Eq. (flU)). Hence, we expect that Eq. alone 
will not be sufficient to reproduce Snla data efficiently, 
at least at low redshift. 



We can do a similar exercise for an evolving effective equation of 
state w B ff(z) instead of constant w c ff. However the corresponding 
expression in this case is much more complicated, so we will not 
present it here. 



III. THE Snla DATA ANALYSIS 

In this section wc shall present a method to place ob- 
servational constraints on the IR corrections ((5j)-([6]) from 
Snla data. We will use the Snla dataset of Hicken et 
at. [13] consisting in total of 397 Snla out of which 100 
come from the new Cf A3 sample and the rest from Kowal- 
ski et al. (l8j . These observations provide the apparent 
magnitude m t h(z) of the Snla at peak brightness after 
implementing the correction for galactic extinction, the 
K-correction and the light curve width-luminosity correc- 
tion. The resulting apparent magnitude m t h(z) is related 
to the luminosity distance cIl(z) = (1 + z)r(z) through 

m th (z) = M(M,H ) + 51og 10 (d L (z)) , (12) 

where M is the magnitude zero point offset and depends 
on the absolute magnitude M and on the present Hubble 
parameter Ho as 

M = M + 5 log 10 (^^j + 25 = M - 5 \og w h + 42.38. 

(13) 

Here the absolute magnitude M is assumed to be con- 
stant after the above mentioned corrections have been 
implemented in mth(z). 

The Snla datapoints are given, after the corrections 
have been implemented, in terms of the distance modulus 



fJ-obs(zi) = m ohs (zi) - M . 



(14) 



The theoretical model parameters are determined by 
minimizing the quantity 



2 _ [MobsQgi) - Mth(z»)] S 

XSnla — 2-^1 



i=l 



a 



(15) 



where N = 397, and c 2 i are the errors due to flux uncer- 
tainties, intrinsic dispersion of Snla absolute magnitude 
and peculiar velocity dispersion. These errors are as- 
sumed to be Gaussian and uncorrclatcd. The theoretical 
distance modulus is defined as 

lM.h(z%) = wth(zi) - M = 51og 10 (di(z)) + , (16) 



where /io = 42.38 — 51og 10 /i and ^i bs is given by Eq. (1141) . 
The steps we have followed for the minimization of 
Eq. (fT5]) in terms of its parameters are described in detail 
in Refs. [TJ-l2l|. 

We will also use the two information criteria known as 
AIC (Akaikc Information Criterion) and BIC (Bayesian 
Information Criterion) , see Ref. [22| and references there 
in. The AIC is defined as 



AIC 



-21n£ n 



2k. 



(17) 



where the likelihood is defined as C cx e x / 2 , the term 
— 21n£ max corresponds to the minimum x 2 and k is the 
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number of parameters of the model. The BIC is defined 
similarly as 



BIC 



-2 In A, 



k In TV , 



(18) 



where N is the number of datapoints in the set under 
consideration. According to these criteria a model with 
the smaller AIC/BIC is considered to be the best and 
specifically, for the BIC a difference of 2 is considered 
as positive evidence, while 6 or more as strong evidence 
in favor of the model with the smaller value. Similarly, 
for the AIC a difference in the range — 2 means that 
the two models have about the same support from the 
data as the best one, for a difference in the range 2 — 4 
this support is considerably less for the model with the 
larger AIC, while for a difference > 10 the model with 
the larger AIC practically irrelevant [22l . [23j . 



IV. RESULTS 

We solve Eq. ©-([5]) numerically to find r(z) for the 
matter-dominated model (SCDM). Then the model is 
tested against the Snla data by using Eqs. (O, ([T5"j). and 
(|16p . Since the parameter 7 does not appear up to third 
order in the expansion of Eq. (|10[) . we will also consider 
the case where 7 = 0. 

In Fig. [T] (left) we present the best fit distance mod- 
ulus versus the redshift z for the best fit ACDM model 
(dashed line), with the present matter density parameter 



"On 



0.289 



+0.023 
-0.022 



) and the SCDM model + IR correc- 



tion of Eq. (fTT7|) with all 3 parameters (solid black line) 
and the 2 parameter case with 7 = (dotted line). For 
the two parameter case we find that a = 1.05^^4 and 
f3 = — 0.51±g J| for ax 2 = 464.031 or a x 2 per degree of 
freedom ~ 1.17, whereas the best fit ACDM has a x 2 per 
degree of freedom ~ 1.18. 

In Fig. [1] (right) we plot the residuals of our best 
fit relative to the ACDM model, i.e. the difference of 
Mbestfit(^) — Macdm(z), for the model with the IR correc- 
tion of Eq. ([5]) with all 3 parameters (solid black line) and 
the 2 parameter case and 7 = (dotted line) . In Fig. [5] 
we show the la and 2a contours for the parameters a 
and (3 of Eq. (|5|) with 7 = 0. The black dot indicates the 
best-fit. 

We should note that the two parameter model fits very 
well the data even if the exact numbers a = 1, /3 = 
— 1/2 and 7 = are used. In this case we find that 
X 2 = 465.462 and a x 2 P er degree of freedom ~ 1.18, 
which is the same as the best fit ACDM model. The 
reason for this success of the model can be seen from 
the first of Eq. (jlip . The standard cosmological model 
predicts values of w c s ~ —1 and Qf^ ~ 0.7, which gives 
a ~ 1. The value of j3 ~ —0.7 derived from the second of 
Eqs. (fn"j) is slightly different from its best fit j3 ~ —0.5 
because of the limitation of the expansion f|T0|) valid only 
in the low redshift regime. 

The original model of Ref. [lU, i.e. Eq. ©, corre- 
sponds to a — /3 — 0, 7 = 1/4. In this case, however, 



the agreement is not very good as the x 2 P er degree of 
freedom is ~ 2.21 with the difference from the ACDM 
model being about 20cr. This discrepancy can be ex- 
plained by the fact that the IR correction does not kick 
in early enough to allow for good compatibility with the 
data, while in low rcdshifts the SCDM behavior of the 
model dominates. This latter property comes from the 
fact that the 7-dependent term appears only at the order 
of z 4 in Eq. ([TO)) . 

If we consider all three parameters a, /3 and 7 to be 
free, then the best fit parameters are a = 0.61, j3 = 1.29 
and 7 = —1.51 for a % 2 = 459.424 or a % 2 per degree 
of freedom ~ 1.17. As it can be seen in Fig. Q] (solid 
black line), the corresponding luminosity distance shows 
a more significant departure from ACDM at high redshift. 

Finally, it is interesting to consider the case where we 
fix the parameters a and (5 to the exact numbers a = 1, 
/3 = —1/2 and allow 7 to vary. In this case we expect 
that by changing 7 we will be able to improve % 2 and still 
be able to compare with ACDM as there will be only one 
free parameter in the model. The result is x 2 = 465.082 
for 7 = 0.052, which is slightly worse (8x 2 ~ 1) than the 
three parameter case, but it is slightly better (5x 2 ~ 0.5) 
than ACDM ( X 2 = 465.513 for O 0m = 0.289). 

In Table|T]we present the results of the two information 
criteria, the AIC and the BIC. According to the AIC the 
three parameter model is the best, with the ACDM hav- 
ing considerably less support, as the difference between 
the two is slightly larger than 2. On the contrary, using 
the BIC we find that the ACDM is now the best and 
also having positive evidence against the other models. 
In all cases, the two parameter model fairs moderately 
with both the AIC and the BIC. 



V. CONCLUSIONS 

We have considered IR modifications of gravity that do 
not imply the presence of a new mass scale in the theory 
and we have studied their compatibility with the Snla 



Table I: Comparison of the one, two and three parameter 
models to ACDM. Note that the differences for the AIC and 
BIC are in both cases with regard to the model with the mini- 
mum value for the corresponding criterion. For the definition 
of the AIC and BIC, see Eqs. (fTTl) and (Tl8)) . 



Model 



AIC 



AAIC 



BIC 



ABIC 



ACDM 

1 parameter 

2 parameters 



467.513 
467.082 
468.031 



3 parameters 465.424 



2.089 471.497 0.000 

1.658 471.066 -0.431 

2.607 475.999 4.502 

0.000 477.376 5.879 
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Figure 1: Left: The distance modulus for the best fit ACDM model (dashed line) and the SCDM model + IR correction of 
Eq. ([5]) with the 3 parameter case (solid black line) and the 2 parameter case with 7 = (dotted line). Right: The residuals 
relative to ACDM for the model with the IR correction of Eq. ([5} with the 3 parameter case (solid black line) and the 2 
parameter case with 7 = (dotted line). The dashed line corresponds to zero. 
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Figure 2: The la and 2a contours in the (a, /3) plane for the 
2 parameter model with 7 = 0. The black dot (a = 1.05, 
/3 = —0.51) corresponds to the best fit. 



data. Our first result is that the mechanism derived in 
Ref. [IH (see also Appendix EJ, Eq. (J]), is not enough, 
by itself, to describe the observed amount of acceleration. 
The absence of free parameters in equation ([T]) (the model 
in [TBI has one parameter less than ACDM) docs not 
make up for the very poor fit of the data. However, a 
more general structure of corrections (U]) can lead to a 
sensibly larger luminosity distance than in the Einstein 
de Sitter universe. In particular, when 7 = in Eq. ([3])- 
©, we have found that the model fits the data very well 
for the values close to the exact numbers a = 1 and f3 = 
— 1/2. It is interesting to consider such sharp numerical 
values, not because of abstract numerology, but because 
a mechanism analogous to that described in Appendix lAl 
very naturally produces coefficients which are integers or 
simple fractions. 

At present it is not clear if the corrections that one 
finds by applying USEP to a field automatically apply 



to, or are inherited by, other types of fields. The sug- 
gested luminosity distance may eventually turn out to 
be produced by considering other fields 7 than the scalar 
field considered in [l5[ or by means of other theoretical 
insights. We also considered the full three parameters 
model ([5])-([H]), whose best fit considerably improves the 
X 2 and is found to be favored over ACDM by the AIC 
but not by the BIC criterion. 

It should be noted that the parameters that we are 
fitting are not coupling constants, and do not appear in 
a Lagrangian. Rather, they are intended as the terms 
of a series expansion that approximate the new "IR- 
complcted" theory starting from GR. As mentioned in 
the introduction, the proposed deformation is present at 
any time during the cosmological evolution; it affects any 
cosmological observation at high redshift, regardless of 
when such an observation takes place, and therefore ad- 
dresses the "coincidence problem" in the most direct way. 

We note that our model requires a rather low value of 
the Hubble constant, Ho ~ 50km/s/Mpc, compared to 
the constraint from the Hubble Key Project with the 
determination of H = 72 ± 8km/s/Mpc [H|. This 
comes from the fact that the Hubble parameter evolves 
as H(z) ~ H Q (1 + zf/ 2 even for z < 0(1) due to the 
absence of a dark energy component. However, methods 
of the determination of Hq that are largely independent 
of distance scales of the Large Mage llanic Cloud Cepheid 
typical ly g ive low values of Hq [26|. For example, Reese 
et al. [23] showed that Sunyaev-Zel'dovich distances to 
41 clusters provide the constraint Hq = 54^3km/s/Mpc 
in the Einstein de Sitter universe. Thus the information 
of Hq alone is not yet sufficient to rule out our model. 



7 It is interesting, for instance, that a different mechanism, based 
on a Casimir-like vacuum energy p3 |. needs a Veneziano ghost 
in order reproduce a density of the right order of magnitude. 
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It would be interesting to sec the effect of the proposed 
IR correction on the angular diameter distance to the 
last scattering surface and estimate the modification to 
the temperature anisotropies in Cosmic Microwave Back- 
ground (CMB). CMB data will be certainly useful to 
place further constraints on our model and to complete 
the picture at higher rcdshift. 
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Appendix A: USEP and the derivation of Eq. (pQ) 



While referring to [15J for details and motivations, it 
is worth summarizing here the basic steps that lead to 
Eq. ([1]) , based on the "ultra strong equivalence principle" 
(USEP). The USEP is a statement about the bare energy 
momentum tensor of the quantized fields on a general 
background, namely: 

USEP: For each matter field or sector suffi- 
ciently decoupled from all other matter fields, 
there exists a state, the "vacuum" , for which 
the expectation value of the (bare) energy 
momentum tensor reads the same as in flat 
space, regardless of the configuration of the 
gravitational field. 

Our starting point is a free scalar field with the action 



s = \J dW^?(3</> 2 - 



in a spatially flat FLRW metric: 

ds 2 = dt 2 -a 2 {f)Ax 2 



(Al) 



(A2) 



The equations of motions of the field and the 00 compo- 
nent of its energy momentum tensor read, respectively, 

d 2 

(j>{t, x) + 3H<f>(t, x) - —\^4>{t, x) + m 2 4>(f x) = 0, (A3) 
a(t) 2 



n{fx) = \ 



h 2 + \(d J cj)) 2 + m 2 (j) 2 



(A4) 



In order to apply USEP we now review the standard cal- 
culation of the energy momentum VEV in such a space- 
time and compare it to the flat space result. 

Upon standard quantization in the Hcisenberg picture 
the field is expanded in creators and annihilators: 



cj){t, x) 



An 



(A5) 

where n is the comoving momentum label in the FLRW 
space; n is a conserved quantity, related to the proper 
physical momentum p by p = n/a(t). From the canonical 
commutation relations [4>(x), tt(x ')} = i(2ir) 3 5 3 (x — x '), 
the commutation relations among the global Fourier op- 
erators are easily derivable, 



{A n ,Al} = 5 3 (n-n'), 



(A6) 



It is customary to choose An as the operator that al- 
ways annihilates the vacuum. The chosen vacuum state 
is therefore implicitly characterized by the choice of the 
mode functions if> n (t), that, by (|A5|) and (|A3|) . satisfy 



0, 



(A7) 



yj n 2 1 a? 



where ui r , 

The solutions of (|A7[) corresponding to the adiabatic 
vacuum [28|, [29[ can be found by a formal WKB expan- 
sion. After quantization, the energy momentum tensor 
of the field (|A4j) becomes an operator whose expectation 
value on the adiabatic vacuum reads [HI, [2j| 



(*,*)> = 



1 



47r 2 a 3 Jo 



H 2 a . _o. 
L0„ + ——+O(n 3 ) 
In 



An . 

(A8) 



The above should be compared to the flat space result 

1 



(T °(t,x)> flat 



47r 2 a 3 J 



n 2 iu n An . (A9) 



It is known that there is a strict connection between 
the geometric properties of a manifold and the spectrum 
of the differential operators [30j or the algebra of func- 
tions [3l| therein defined; such abstract characterizations 
have occasionally been used for generalizing common ge- 
ometrical concepts and the description of spacetime it- 
self [3l|, However, so far, attempts in this direction 
have always been applied to the UV and intended to mod- 
ify spacetime at the smallest scales. Here we would like 
to modify the IR-spcctral properties of the FLRW metric 
(and therefore its geometry) in order to enforce USEP. 
The proposed deformation is argued to correspond to a 
breakdown of the metric Riemannian structure at dis- 
tances comparable to 

So the idea is to modify the physics in the infra-red but 
strictly maintain the equations and the relations valid lo- 
cally such as the field equations (|A3|) and the form of the 
energy momentum tensor (|A4[) . We choose a point in 
FLRW (x = 0) and make a formal Taylor expansion of 



7 



which GR is the zeroth order. In the spirit of a general 
spectral deformation, we conjecture a mismatch between 
the "metric-manifold" Fourier labels n and the physical 
momenta k that locally define the infinitesimal transla- 
tions and the derivatives of the local fields. In other 
words, we now write the field in x ~ as 

0(t,2«O) = j dn 3 [^(tje^+^ftje-*^] . 

(A10) 



where 



H 2 a 2 , . , 
1 — ^ — h higher order 



(All) 



Note that, when derivatives of the field are taken in 
x = 0, a factor of k, instead of n, drops. The form of 
above relation, which is one of the main results of [l5j . 
is dictated by the request that the quadratically diver- 
gent time dependent piece of (|A8j) disappears. In other 
words, that is the first order correction in order to im- 
pose USEP upon this particular GR solution. The cor- 
rected mode equation is in fact obtained by substituting 
(lAlOp into (|A3|) . which is assumed to be strictly valid be- 
cause it applies locally. To the modified mode equation, 
i> k {t) + 3Hip k (t) + u>lipk(t) = 0, the same WKB expan- 
sion can be applied and the quadratic divergence in (|A8|) 
is reabsorbed just by re-expressing (0| Tq|0) in terms of 
the appropriate "flat measure" time-independent Fourier 
coordinates n: 



(T°(t,xn0)) = 



1 



At: 2 a 3 



H 2 a 
~2k~ 



0(n 



In order to define local quantities away from the origin 
we exploit the assumption of spatial homogeneity and use 
the translation operator 



T(A) = cxp ( -iX ■ / d 3 n k(n) A\A n 



(A12) 



that we obtain by simple exponentiation of the (mod- 
ified!) momentum operator. The parameter A is the 
comoving distance. The field away from x — is thus 
defined as 0(i, A) = T,(A) <j>{t, 0) Tr x {\) and reads 



0(i,A) 



<Tn 



Mty~ k - x A H +r k (t)e-' rk - x Ai 



(A13) 

From the above expression it is straightforward to cal- 
culate the modified commutator between the canonical 
momentum 7r(0) = a 3 0(O) and the field at comoving dis- 
tance A, 



[7r(O),0(A)] = -i U 3 (A) + 



1 H 2 a 2 
8^ A 



(A14) 



Note that there is a potential ambiguity in defining the 
time derivative of a displaced operator. By deriving 
(|A13| we get 



n(X) = a 3 I d 3 n \^ k (t)e ik - x A H + ^Je"**"^ 

+ ia 3 fd 3 n(k-X)- Ufe e ttX A rl -r k e"* x 4 



The second line in the above equation is there because k 
is time dependent. However, if we just apply the transla- 
tion to 7r(0), instead of deriving 0(A), those terms would 
not be there. Therefore, for consistency, we need to 
make them ineffective at the required order of approxi- 
mation. This can be done by imposing that [0(0), 7r(A)] = 
— [7r(0), 0(A)]. Because of the second line the last equa- 
tion, the commutator between 0(0) and 7r(A) gives 



[0(O),7T(A)] 



-[tt(O),0(A)] 
-2i 



(2tt) 3 



d 3 ne" 



ik-X 



\4>n\ 2 (k-\)\ 



(A15) 



the last term being the spurious contribution. In order 
to get rid of it, we have to make the comoving physical 
distance A also time dependent. This effectively means 
that, after an infinitesimal time step di, we have to re- 
consider the field translated, from x ~ 0, not by the same 
comoving distance A, but by a slightly different amount. 



At high momenta/small distances, since \ipk 
dn the integral in the second term of (|A15[) reads 



1/n, 



1 f°° 



d 3 ne^™ X - 



H 2 a 2 
2n 2 



,(H 2 a 2 ) 



2n 2 



We make the ansatz A = b\ 2 (H 2 a 2 )' A, where b is a num- 
ber to be determined. We get 



d d ne-* fc ' A hM (fc-A)' 



i{H 2 a 2 } 



da 



d 3 n 



bX 2 



2n 3 ) 



-ia n-X 



a=l 



The last integral can be regularized by setting n~ 3 — > 
n~ 3+£ and taking the e — > limit only after deriving 
with respect to a. The result is null for b = 1/4, which 
fixes the time dependence of A: 



A = A 



3 (H 2 a 2 ) 



(A16) 



Appendix B: Exact solutions 

Here we present some analytical solutions for the dif- 
ferential equation If 7 — we obtain the following 
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analytic solution: 
H r(z) = 



4/yr 



-1 - 2a + V5 + 2y/S/[(l + z)^/ 2 - 1] ' 

(Bl) 

where 5 = (1 + 2a) 2 - 16/3. 

When 7^0, finding the solution is much more difficult 
and can only be given in an implicit form. For example, 
let us consider the correction in Eq. $2$, i.e. a = /? = 
and 7 = 1/4. Setting R(z) = r(z)y/l + z, we get a 
differential equation for the function R(z): 



dR(z) 
dz 



R(z) 



R{zf 



1 + z 2(1 + z) 4(1 + z)' 



(B2) 



with initial conditions R(Q) = and (di?/dz)(0) = 1. By 
direct differentiation it can be shown that the solution to 



Eq. (|B2I) is given in an implicit form by: 



R(z) 
x 3 



V ^2 



(B3) 



where the parameters X\,X2, x 3 are the roots of the poly- 
nomial equation: 

4 + 2 Xl + xf = . (B4) 

When all three parameters a, /3, and 7 are not zero, 
then we can still find an implicit solution for r(z) but in 
this case it is very complicated, so we will not reproduce 
it here. 
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